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Abstract — We analyze the capacity scaling laws of clustered 
ad hoc networks in which nodes are distributed according to 
a doubly stochastic shot-noise Cox process. We identify five 
different operational regimes, and for each regime we devise 
a communication strategy that allows to achieve a throughput to 
within a poly-logarithmic factor (in the number of nodes) of the 
maximum theoretical capacity. 

I. Introduction and related work 

The capacity of ad hoc wireless networks has been 
traditionally studied considering single-user communication 
schemes over point-to-point links [1|. Only recently Q, Q, 
(4), information-theoretic scaling laws of ad hoc networks 
have been investigated, showing that multi-user cooperative 
schemes can achieve much better performance than traditional 
single-user schemes, especially in the low power attenuation 
regime. 

In this paper, we follow the stream of work 0, 0, 
J4), analyzing the information-theoretic capacity of clustered 
random networks containing significant inhomogeneities in 
the node spatial distribution. In particular, we consider nodes 
distributed according to a doubly stochastic Shot-Noise Cox 
Process (SNCP) over a square region whose edge size can 
scale with the number of nodes. 

We provide both information-theoretic upper-bounds to the 
achievable capacity and constructive lower-bounds, which are 
asymptotically tight to within a poly-log factor (in the number 
of nodes). Our study reveals the emergence of five operational 
regimes, in which different communication schemes combined 
with proper scheduling/routing strategies must be adopted to 
achieve the system capacity. 

With respect to previous work, we provide several contri- 
butions. First, the analysis in 0, Q is limited to networks 
in which nodes are uniformly distributed. In contrast, our 
complete characterization of the network capacity achievable 
under the SNCP model extends the analysis to a much broader 
class of network topologies (including the uniform distribution 
as a special case), which can take into account the clustering 
behavior usually found in real systems. 

Second, the impact of inhomogeneities in the node spatial 
distribution has been first investigated in |4), where authors 
have found that for small path-loss exponents (i.e., a £ (2, 3]) 
the capacity does not depend on how nodes are placed over 
the area. Instead, they show that capacity is significantly 
affected by the network topology for large path-loss exponents 



(i.e., a > 3). However, their characterization of the capacity 
achievable for large path-loss exponents is limited to the 
case of adversarial node placement under a deterministic 
(given) degree of network regularity. Moreover, they impose a 
minimum separation constraint between the nodes which does 
not allow to introduce highly dense clusters over the area. 
At last, the analysis in [4] is limited to the case of extended 
networks (i.e., networks whose area grows linearly with the 
number of node£3). 

Third, our constructive lower bounds require to employ 
novel scheduling/routing strategies in combination to existing 
cooperative communication schemes. Such strategies represent 
an important contribution in themselves, as they could be 
adopted to cope with the nodes spatial inhomogeneity in more 
general topologies which cannot be described by the SNCP 
model considered here. 

At last we emphasize that this work extends 0, (6), where 
we have analyzed the capacity of networks in which nodes are 
distributed according to a SNCP model, but considering single- 
user communication schemes only (i.e., traditional point-to- 
point links). 

II. System Assumptions and Notation 
A. Network Topology 

We consider a network composed of a random number N 
of nodes (being E[N] = n) distributed over a square region 
O of edge length L, where L takes units of distance. The 
network physical extension L scales with the average number 
of nodes, since this is expected to occur in many growing 
systems. Throughout this work we will assume that L = rO ', 
with 7 > 0. To avoid border effects, we consider wrap-around 
conditions at the network edges (i.e., the network area is 
assumed to be the surface of a bi-dimensional Torus). 

The clustering behavior of large scale systems is taken into 
account assuming that nodes are placed according to a shot- 
noise Cox process (SNCP). An SNCP Q over an area O 
can be conveniently described by the following construction. 
We first specify a homogeneous Poisson point process C of 
cluster centres, whose positions are denoted by C = {cj}jL lt 
where M is a random number with average i£[M] = m. Each 
centre point Cj in turn independently generates a point process 

In (3J authors recognized the importance of letting the network area scale 
with the number of nodes in a general way, as this gives rise to a richer set 
of operational regimes. 



of nodes whose intensity at £ is given by qk(cj,£), where 
q G (0, oo) and k(cj, •) = c j\\) is a rotationally invariant 

dispersion density function, also called kernel, or shot; i.e., 
k(cj,-) depends only on the euclidean distance ||£ — Cj\\ of 
point £ from the cluster centre Cj. 

Moreover we assume that fc(||£— Cj || ) is a non-negative, non- 
increasing, bounded and continuous function, whose integral 
J a k(cj,£) d£ over the entire network area is finite and equal 
to 1. In practice, the kernels considered in our work can 
be specified by first defining a non-negative, non-increasing 
continuous function s(p) such that p s(p) dp < oo and 
then normalizing it over the network area O: 
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Notice that, in order to have finite integral over increasing 
network areas, functions s(p) must be o(p~ 2 ), i.e., they 
must have a tail that decays with the distance faster than 
quadratically. In the following, we will be especially interested 
in functions s(p) whose tail decays as a power-law: 



s(p) = min(l,p-*) for S > 2, 



although our results apply to more general shapes as well. 
At last, in our asymptotic analysis we can neglect the 



normalizing factor J„ s(||£ - 



idc = e(i) 



Under the above assumptions on the kernel shape, quantity q 
equals the average number of nodes generated by each cluster 
centre (all cluster centres generate on average the same number 
of nodes). In our work, we let q scale with n as well (clusters 
are expected to grow in size as the number of nodes increases). 
This is achieved assuming that the average number of cluster 
centres scales as m = n v , with v G (0, 1). Consequently, the 
number of nodes per cluster scales as q = n 1 ^". 

The overall node process Af is then given by the super- 
position of the individual processes generated by the cluster 
centres. The local intensity at £ G O of the resulting SNCP is 

M 
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Notice that $(£) is a random field, in the sense that, 
conditionally over all (cj), the node process Af is an (inho- 
mogeneous) Poisson point process with intensity function $. 
We denote by X = {Xi]f =l the collection of nodes positions 
in a given realization of the SNCP. 

Let d c — L/y/m = n~ l ~ 1 '/ 2 be the typical distance between 
cluster centres. More precisely, d c is the edge of the square 
where the expected number of cluster centres falling in it 
equals 1. We call 

cluster-dense condition the case 7 < v/2, in which d c 

tends to zero an n increases; 

cluster-sparse condition the case 7 > v/2, in which d c 

tends to infinity an n increases. 
Figure [TJ shows two examples of topologies generated by 
our SNCP, in the case of n = 10,000 and 7 = 0.25. In both 
cases we have assumed s(p) = min(l, p~ 2 ' 5 )- The topology in 
Figure [T(a)1 has been obtained with v = 0.6, hence it satisfies 
the cluster-dense condition (7 < v/2). The topology in Figure 
I l(b)| corresponds to v = 0.3, and provides an example of the 
cluster-sparse condition (7 > v/2). 




(a) SNCP with v = 0.6. 



(b) SNCP with v = 0.3. 



Fig. 1. Examples of topologies with n = 10, 000 nodes distributed over the 
square 10 X 10 (7 = 0.25). In both cases s(p) ~ p -2 ' 5 . 



Recall that the local intensity of nodes at point £ can be 
written as <£>(£) = ^ . q k(cj , £). We define the quantities: $ = 
su P£eo an d = infseo The following lemma, 

proven in [5|, characterizes the asymptotic behavior of $ and 
(1) !: 

Lemma 1: Consider nodes distributed according to the 



SNCP. Let rj{m) = d c y/\og m. If 77(771) = o(l), it is possible 
to find two positive constants gi,Gi with g\ < G\ such that 

V£ g O 

<m< Gi^ w.h.pH (2) 

When 77(771) = 0(1), it is possible to find two positive con- 
stants <?2,G2, such that, w.h.p., $ > <72<? log m s(d c y / log m) 
and $ < G^q logm. 

The above result implies that $ = 0($) in the cluster- 
dense condition, i.e., when 7 < v/2 (which implies d c = 
o(l/y / logm)), whereas $ = o($) in the cluster-sparse 
condition, i.e., when 7 > v/2 (which implies d c = ui(l)). 

B. Communication Model 

We use the same channel model as in 0, 0, 0. Consider 
the generic time t, and let V(t) be the set of nodes transmitting 
at time t. The signal received at time t by a node k is 

y k [t] = Kk[t]xi[t]+z k [t] 
iev(t)\{k] 

where Xi[t] is the signal emitted by node i, and {z/c 
are white circularly symmetric Gaussian noise, independently 
and identically distributed (i.i.d.) with distribution A/c(0, iVo) 
(with zero mean and variance No per symbol). The complex 
baseband-equivalent channel gain between i and k at 

time t is 

hi, k [t] = VGdl 

where G is a constant gain, a > 2 is the path-loss exponent, 
and {f9ife[rj]}i.fc are i.i.d. random phases with uniform distri- 
bution in [0, 2tt), which are assumed to vary in a stationary 
ergodic manner over time (fast fading). Moreover, 
and {dik}i,h are a ls° assumed to be independent, Vi, k. We 
should mention that a recent work J8 | has put in discussion 
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2 Throughout the paper, we adopt the terminology 'with high proba- 
bility' (w.h.p.) to indicate events/properties that occur with a probability 

P ' 
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TABLE I 

System parameters (n.a. = not applicable) 



Symbol 


Definition 


scaling exponent 


L 


edge length of network area 


7 > 


m 


average number of clusters 


< v < 1 


P 


per-node power budget 





a 


path-loss exponent 


n.a. 


8 


decay exponent of s(p) 


n.a. 


d c 


typical distance between cluster centres 


7-1//2 


q 


average number of nodes per cluster 


1-u 



the validity of this assumption for very large n in the low 
path-loss regime a £ (2, 3). However, the strong impact of the 
assumptions on the location of scatterers suggests that channel 
modelling in the low path-loss regime for very large networks 
is somewhat delicate and requires further investigation. 

We assume that each node is source and destination of a 
single flow, and that the resulting N flows (with E[N] = n) 
are established at random without any consideration of node 
locations. Let A(n) be the largest uniformly achievable rate 
of communications between sources and destinations. The 
aggregate system capacity is C(n) — n\(n). At last, we 
impose an average power constraint of P on the transmissions 
performed by each node, where P is a constant. 

Table Q] summarizes the parameters of our model. For the 
quantities that are allowed to scale with n we have reported, in 
the third column, the restrictions on the scaling exponent in n, 
i.e., the assumptions on log„(<parameter>). Note that d c and 
q are not native parameters, since they are derived from others, 
however we have included them in the table for convenience. 

III. Summary of Results 

Similarly to previous work Q, 0, we express our results 
in terms of the scaling exponent ec of the network capacity, 
defined as, 

logC(n) 
e c = km — 

rn-oo log n 

The scaling exponent allows to ignore all poly-logarithmic 
factors, i.e., factors which are 0(logn) fe , for any finite k. 
Since our lower and upper bounds differ at most by a poly- 
logarithmic factor, the corresponding scaling exponents match, 
and we can claim that our characterization of the network 
capacity in terms of the scaling exponent is exact. 

Results are reported in Table [XT] The scaling exponent 
takes different expressions as functions of the four system 
parameters {a, 7, 6, v}, under the conditions specified in the 
third column of Table HTl In particular, we can distinguish five 
operational regimes, denoted by latin numbers 1,11,..., V, as 
reported in the second column of the table0. It can be verified 
that ec varies with continuity in the four-dimensional space 
of parameters {a, 7, 6, v). We observe that, under any regime, 
ec is a non-increasing function of parameters {a, 7, 5} and a 
non-decreasing function of v. Figure [2] provides a graphical 
representation of the results in Table [ID for the particular case 
of v = 0.3 and 8 = 2.5. 

3 In the last two rows of the table, the actual regime depends on which term 
prevails in the max(-) expression used in column one: we are in regime III 
if the first term is bigger than the second one. 




Fig. 2. Scaling exponent of network capacity as function of a and 7, for 
v = 0.3 and 8 = 2.5. Different marks are associated to the five possible 
regimes. 

IV. Upper Bounds 

Upper-bounds are obtained extending the approach in 0, 
0, El, which is based on the computation of a bound to 
the information flow passing through a cut that divides the 
network in two parts. 

First, by leveraging percolative arguments (see 0), it 
is possible to find a strip of width A, with A = 
Q((qs(d c ) logn) -1 / 2 ) in the cluster-sparse condition (i.e., for 
7 > v/2) and A = Q(L/y/n) in the cluster-dense condition 
(i.e., for 7 < v/2), which divides the network area in two 
parts, and satisfies the following properties: i) the considered 
strip is empty of nodes; ii) every cluster centre lies at a distance 
greater than g d c from the strip, for a sufficiently small constant 
9- 

Then, the information flow C(S,T>) from sources S on the 
left to destinations T) on the right can be bounded by the power 
transfer Vs,t> through the strip, as in 0, 0, El- Let dn, be 
the euclidean distance between node i and node k. According 
to El: 

C(S,V) < bV s , v = b ]T PdT k a (3) 
ies,keT> 

for any / x 
b > 4 max 1 , max > — — 

being, for every k, ^ ie5 ( ''o ) = 0(k>g 5 n). 

To estimate Vs,v, the left and right domains are partitioned, 
respectively, into squarelets {Afcjfc and {Bh]h, obtaining: 

V s ,v = P J2 d ^ < p J2J2^kU(A k )U(B h ) 

i£S,keV h k 

where d hk is the minimum distance between points of A k 
and points of B>h, while function U(Ak) (U(Bh)) provides an 
upper bound to the number of nodes in A k (B^). To obtain 
tight upper bounds the size of A k and Bh must be carefully 
chosen since, by increasing the size of A k and Bh, on the one 
hand we obtain tighter bounds for U(Ak) and U(Bh)', on the 



TABLE II 

Scaling exponent of network capacity = 1 — v — 5(7 —v/2). 





regime 


conditions 


1 

2 — c*7 
Q-I-Q7 

a-2 

max [2 - a 7 + (a - 3)| , 7 + /3§Ei] 
max [2 - q 7 + (a - 3) | , 7 + 


I 
I 

II 

III or IV 
III or V 


cry < 1 

«7 > 1 A a < 3 

«7>1Aq>3A > 7 - f 
«7>lAa>3A <7-fA/3>0 
«7>lAa>3A ^2 <7-fA/3<0 



other, we obtain looser bounds d hk for the distance between 
nodes in A k and nodes in Bh- 

Furthermore, when A = o(l), a tighter bound can be 
obtained applying the Hadamard inequality (see [3]) to extract 
from the information flow the contribution of destinations 
receiving signals whose strength diverges. This contribution, 
associated to nodes in proximity of the cut, is in turn evaluated 
applying the Hadamard inequality iteratively, so as to split 
it into the contributions associated to individual destinations 
(which can be interpreted as MISO systems running in par- 
allel). Each individual contribution is then bounded applying 
similar arguments as in |3l. 

The above mentioned five regimes derive from the fact that 
the dominant contribution to C(S,T>) changes while varying 
the system parameters. In regime I the dominant contribution 
is due to nodes lying at distance 0(X) from the cut; in regime 
II the dominant contribution is provided by nodes which are 
jointly at distance uj(d c y/logn) and o(L); in regime III it is 
due to nodes at distance 0(d c ^/\ogn) and f2(d c ); in regime 
IV it is due to nodes at distance o(d c ) and uj(1/t/W); at last, 
in regime V the dominant contribution is provided by nodes at 
distance Q(1/t/W). A detailed derivation of the upper-bounds 
can be found in j9j. 



V. Lower Bounds 

For each operational regime, it is possible to devise a 
communication scheme that approaches the corresponding 
upper bound to within a poly-log factor. All of our proposed 
schemes work as follows: first, a subset of nodes is identified, 
which forms the main infrastructure through which data is 
transferred across the network area. A finite fraction of time 
is then assigned to the rest of the nodes to exchange traffic 
with the nodes belonging to the main infrastructure (if needed). 
More precisely, time is divided into regular frames, each one 
comprising three phases of equal duration: i) an access phase, 
in which sources not belonging to the main infrastructure send 
data to the infrastructure; ii) a transport phase, in which data 
is transferred over the infrastructure; iii) a delivery phase, 
in which data is sent from the infrastructure to destinations 
not belonging to it. Since the delivery phase is analogous 
to the access phase (by exchanging the role of transmitters 
and receivers), we will focus on the access phase only, after 
presenting the transport phase. 

Before proceeding, we report the lower bounds obtained in 
O for homogeneous networks. Given a Homogeneous Poisson 
Process (HPP) of intensity ip over a square (or disc) of edge 
(radius) L, it is possible to achieve the aggregate capacity 
C„(L, ip, a): 



u}(N 2 - e L- a ) 



N >L a 
N < L a , a < 3 



u}(N- e Li/j—) N < L a , a > 3, V = w(l) 



u(N- e LiP^) 



a > 3, V = 0(1) 



(4) 

w.h.p. for any e > 0. In the above expressions N = %/jL 2 is 
the average number of nodes in the system. 

A. Transport phase 



concerns the transport phase, our proposed 
can be considered as special cases of a general 



For what 
schemes 

class of scheduling-routing strategies, according to which 
the network area is partitioned into cells of edge size A 
cooperative multi-hop strategy is applied, in which MIMO 
communications are established between the nodes belonging 
to neighboring cells, and global multi-hopping at the cell 
level is employed to transfer data through the network. The 
proposed schemes essentially differ in: i) the subset of nodes 
which are used as the main infrastructure; ii) the chosen value 
of the cell edge size I. In particular, the value of I allows us to 
classify our schemes into five main communication strategies 
(for the transport phase) which can be associated by a one-to- 
one correspondence to the five operational regimes reported in 
Table HIl 

I: global MIMO, in which I = Q(L), and nodes employ 

a MIMO communication scheme at global network scale, 

without the need of cell multi-hopping; 

II: cooperative super-cluster hopping, in which nodes 

employ a cooperative multi-hop scheme, where I = 

w(d c Vlogn); 

III: cooperative inter-cluster hopping, in which I = 
<d(d c \/\ogn), i.e., the cell edge size is closely related to 
the typical distance d c between cluster centres; 
IV: cooperative sub-cluster hopping, in which I = o(d c ) 
and I = w(l/$), i.e., the cell edge size is smaller (in 
order sense) than the typical distance between cluster 
centres, yet the cell is large enough to allow cooperation 
among an increasingly number of nodes falling in it; 
V: traditional multi-hop scheme, in which I = 
0(l/v / $), and nodes resort to the traditional point-to- 
point multi-hop scheme, since there is no advantage (in 
order sense) in employing cooperative techniques. 
Notice that the above five strategies for the transport phase 
are applied to different infrastructures, which are selected 
depending on the combination of system parameters. The basic 
tool that we use to extract a subset of nodes forming the main 
infrastructure is a standard thinning technique, that can be 
applied to our class of point processes in the sense specified 
by the following lemma. 



Lemma 2: Consider nodes X = {^}^ placed according 
to the considered SNCP. Then a subset of nodes ZCX can 
be found w.h.p. such that Z forms a homogeneous Poisson 
process with intensity J> , where $ = gify in the cluster- 
dense condition and $ = g? q log m s (d c ^/log m) in the 
cluster-sparse condition. Here g\ and <?2 are the constants 
defined in Lemma Q] 

We identify the following three main infrastructures: 

dense infrastructure, which is used in regimes I and II, 
but only for the cluster-dense condition (7 < v/2). In 
this case, we can apply Lemma [2] and extract a subset Z 
of cardinality 9(n), which can sustain the same capacity 
of a homogeneous system with n nodes; 
clusters-core infrastructure, which is used in regimes 
I, II, III, for the cluster-sparse condition (7 > v/2). In 
this case, the set Z is formed by all nodes falling within 
a finite distance from their cluster centre. The cardinality 
of this set is still 0(n); 

sparse infrastructure, which is used in regimes IV and 
V, for the cluster-sparse condition (7 > v/2). In this case, 
we can apply Lemma [2] and extract a subset Z of points 
with density $ = 6(n' 3 ), where (3 = l-v-5(-)-v/2). 
The cardinality of this set is o(n). 
Since both the dense infrastructure and the sparse infrastruc- 
ture form a HPP, their capacity can be immediately obtained 
applying existing results for homogeneous system. The cluster- 
core infrastructure is not a HPP, however it can be regarded 
as being uniformly dense at resolution higher than d c . Since 
in regime I,II,III the cell edge size is Vl{d c yJ\og n), MIMO 
communications between cells occur as if nodes in Z were 
uniformly distributed (see |9| for more details). Moreover, it 
can be shown that the clusters-core infrastructure can sustain 
the load due to the cooperation overhead required within each 
cell, but we omit the details here. 

B. Access phase 

We recall that the access phase is used by sources to inject 
their traffic over the main infrastructure. Since the system 
capacity is ultimately determined by the main infrastructure, 
the goal is to design an access phase that does not constitute 
a system bottleneck, while at the same time inducing a 
uniform traffic matrix over the main infrastructure. These 
design principles led us to select the following three access 
strategies: 

SISO access scheme. This is the simplest strategy, and 
it is used to access the dense infrastructure. In this 
case, it is sufficient to employ a single-hop point-to-point 
transmission between each source and one of the closest 
nodes belonging to Z, thanks to the fact that the network 
is almost uniformly dense; 

SIMO access scheme. This is used to access the nodes 
of the clusters-core infrastructure, employing a SIMO 
technique similar to the relaying scheme proposed in |40; 
hierarchical access scheme. This is used to access the 
nodes of the sparse infrastructure, and required us to 

4 In (4j, authors present a technique that allows nodes located in low- 
density areas to relay their data over densely populated areas, by exploiting 
the diversity gain intrinsically available in high-density regions thanks to the 
presence of many nodes acting as an array of receiving antennas. 



develop a novel scheduling-routing strategy specifically 

tailored to this case. 
Due to lack of space, we restrict ourselves to a brief 
description of the hierarchical access scheme, which is the 
most intriguing on^fl In this case, traffic produced within 
highly dense regions of the network area (e.g., the clusters 
cores in Figure [T(b)| needs to be gradually spread out through 
a sequence of intermediate, local transport infrastructures 
nested one within the other, This construction is needed both 
to avoid the formation of local bottlenecks around the cluster 
centres, and to evenly balance the traffic towards the node of 
the main infrastructure. Intermediate transport infrastructures 
are obtained by applying the thinning technique of Lemma 
12 within certain domains (specified later), surrounding the 
clusters' centres, nested one within the other. To simply and 
effectively balance the traffic data are delivered within each 
local infrastructure to randomly destination nodes. 

The sequence k = 0, 1, . . . , K max of nested domains is 
carefully chosen in such a way that: i) the first domain in 
the sequence coincides with the network area, hence the cor- 
responding infrastructure is the main transport infrastructure 
of the network, of density $, which is shared by all data flows; 
ii) the infrastructure extracted in each domain k > can pass 
to the infrastructure of domain k — 1 all traffic generated by 
nodes contained in it; iii) the total number of domains grows 
at most like log n. 
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Fig. 3. Example of construction of nested domains Ok for the topology 
depicted in Figure [1(b)] Domain 0\ is characterized by d± = 0.5d c . 



Conditions i) and ii) guarantee that the system capacity 
is throttled by the lowest infrastructure (the main transport 
infrastructure) and no bottleneck arises within any higher 
infrastructure. Condition iii) guarantees that, even if we devote 
to each layer-fc infrastructure the same fraction of time, the 
total overhead due to the access phase causes at most a log n 
loss in the overall system capacity. 

We now specify one possible way to jointly achieve the three 
conditions above. We build a sequence of nested domains Ok, 
k = 0, 1, 2 ... , i\~ m ax, as follows. The first domain is Oq = O, 
meeting condition i). 

For the generic point £ 6 O, let d mm (£) = mirij c,-|| be 
the distance between £ and the closest cluster centre. We define 
domains Ok, for k > 1, as follows: Ok — {£ € O : d m in(£) < 
dk}, where dk are a set of decreasing distances, i.e., di > c?2 > 
. . . > dK max ■ Domain Ok is, in general, composed of a random 

5 The interested reader is referred to (5) for a detailed description and 
analysis of all access schemes. 



number Jk of disjoint regions (Jk < M), corresponding to 
the connected components of the standard Gilbert's model of 
continuum percolation iflOl with ball radius dk- Figure[3]shows 
examples of domains Ok having different values of dk- Let 
be the set of disjoint regions (1 < j < Jk) forming 
domain Ok- 

We set the largest dk, namely d\, equal to d\ — fid c , 
where [i is a small constant. Choosing [i sufficiently small, 
in such a way that the associated Gilbert's model is below 
the percolation threshold (we need /i < \i* , where \i* s» 0.6), 
we have the property that the maximum number of clusters 
centres belonging to the same region If is O(logn) w.h.p. 
[ 10 1 . Since by construction Ok+i C Ok, the same property 
holds for all k > 1. It follows that, in terms of physical 
extension, the area \Tl\ of region 2| lies w.h.p. in the interval 

nd 2 k < <7rd|logn. 

We further observe that the density of nodes at any point 
within Ok (k > 1) can be lower bounded by A^ = q d k s , by 
considering the contribution of the closest cluster centre only. 
Hence, it is possible to extract from Ok (fc > 1) a set of points 
Zfc forming a HPP with intensity A^. Note that in the domain 
Oq we have Ao = $. Distances dk, for k > 2, are then as- 
signed in such a way that A& = 2 1 Ai, i.e., the intensities of 
the nested transport infrastructures form a geometric progres- 
sion. This requires to set dk = d\l ~ . Since the maximum 
node density in the network is $ < GiQ logm (see Lemma 
CO, we have K max = 1 + |_log 2 (<3Togm/Ai)J = O(logn), 
hence the total number of domains satisfies condition iii). 

It remains to show that each domain k < K max can receive 
the traffic generated by domain k+1. To this purpose, we need 
to show that each region II can handle the traffic produced 
by all components of domain k + 1 nested in it. Let T-L k +i ^ e 
the set of indexes h of regions T^+x filing m 2f • Moreover, 
let Ml be the number of cluster centres falling within I?. 

The area of II can be expressed as \IV\ — M J k nd\C,k, 
where (k < 1 is a reduction factor that accounts for the 
overlapping among the discs of radius dk forming region I?. 
The sum of the arecis of till nested regions is instead 

given by £ fc£K i l Z fe+il = M l^d 2 k +iCk+i, where Cfc+i > Cfc 
because the degree of overlapping among the discs reduces 
for decreasing values of dk- Since (dk/dk+i) 2 — 2 2 / 5 , we 
conclude that the ratio between \I 3 k \ and J2heH j l-^fe+ll * s 
bounded. This is important, as it allows to exploit to full 
capacity of the infrastructure extracted in Ti to spread out 
the traffic coming from nested regions l k+1 over the larger 
region 2|. 

Moreover, using the expressions it can be shown that 
the aggregate capacity of nested regions X k , l is larger than 

the capacity of region II. This allows to conclude that domain 
Oo (i.e., the main infrastructure) acts as the system bottleneck. 
Indeed, the number of points in li is 

MlnXkdlCk = M^7rA 1 d 1 2( fe - 1 )( 1 -f)C fe 

The total number of points in regions l^+i rias me same 
expression, substituting k with k + 1. Since S > 2, and 



Cfc+i > Ck, the total number of points in regions is 
larger. This guarantees that the aggregate capacities of the 
nested infrastructures is higher than the capacity of II in the 
first regime of in which C n = uj(N 1 ~ e ). 

In the third regime of (0J, the capacity (either of region 
I 3 k or the aggregate capacity of nested regions l(f +1 ) would 
be proportional to 2 fe[ ^~7- £ ( 1 -f)l£ fc . Since > 1 > 

i, and e is small, the capacity increases with k. At last, in 
the forth regime of (0]l the capacity would be proportional 
to 2 fc [T--^(i-§)]^, which again increases with k. One 
can verify that capacities still form a non-decreasing sequence 
when we change regime passing from layer k to layer k + 1. 

We conclude that the chosen sequence of nested local 
infrastructures satisfies the conditions that allow to balance 
the traffic towards the nodes of the main infrastructure at most 
with a log n penalty factor to the system capacity. 

VI. CONCLUSIONS 

We have characterized the asymptotic capacity of networks 
whose nodes are distributed according to a doubly stochas- 
tic shot-noise Cox process. This point process provides an 
interesting, analytically tractable model of clustered random 
networks containing large inhomogeneities in the node density. 
Our study has revealed the existence of additional operational 
regimes with respect to those identified in previous work, 
and the need of novel scheduling and routing strategies, 
specifically tailored to each regime, to approach the maximum 
system capacity. 
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